Neutron and Proton Transverse Emission Ratio Measurements and the Density 
Dependence of the Asymmetry Term of the Nuclear Equation of State 
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Recent measurements 

112,124c„ i 112.124 



of pre-equilibrium neutron and proton transverse emission from 
Sn+ 1-La,J-M Sn reactions at 50 MeV/A have been completed at the National Superconduct- 
ing Cyclotron Laboratory. Free nucleon transverse emission ratios are compared to those of A=3 
mirror nuclei. Comparisons are made to BUU transport calculations and conclusions concerning the 
density dependence of the asymmetry term of the nuclear equation-of-state at sub-nuclear densities 
are made. The double-ratio of neutron-proton ratios between two reactions is employed as a means 
of reducing first-order Coulomb effects and detector efficiency effects. Comparison to BUU model 
predictions indicate a density dependence of the asymmetry energy that is closer to a form in which 
the asymmety energy increases as the square root of the density for the density region studied. A 
coalescent-invariant analysis is introduced as a means of reducing suggested difficulties with cluster 
emission in total nucleon emission. Future experimentation is presented. 



PACS numbers: 25.70.-z,25.70.Mn,21.65.+f,26.62.+c 



The nuclear symmetry energy increases the masses of 
nuclei with very different neutron and proton concentra- 
tions and limits the neutron concentration and maximum 
neutron number N of any element. In the interior regions 
of neutron stars, where neutrons may comprise over 90% 
of the matter, the symmetry energy may contribute the 
bulk of the pressure supporting the star |]J . The nuclear 
equation-of-state at densities of 0.5< — <10 (where po 
is the nuclear saturation density) governs many of the 
neutron star macroscopic properties, including radius, 
moment of inertia, core structure 0], cooling rates, and 
thepossible collapse of a neutron star into a black hole 

□ DQ 

Constraints on the symmetry energy at sub-saturation 
density have been obtained from measurements of the 
diffusion of neutrons and protons between nuclei of dif- 
ferent asymmetry S = (N-Z)/(N+Z) in peripheral colli- 
sions 

\Hm ( Here N and z 

arc the relevant neutron and 
proton number of the nuclei.) Measurements of nuclear 
masses, isovector collective excitations and neutron skin 
measurements may also provide constraints at densities 
less than pq. Nevertheless, the density dependence of the 
symmetry energy is not known well enough to constrain 
the relevant neutron star properties. Using an alterna- 
tive approach to improve present constraints, new mea- 
surements of the ratios of neutron and proton spectra in 
central heavy ion collisions are presented and compared 
to transport theory calculations. These calculations dis- 
play a strong sensitivity to the density dependence of 
the symmetry energy, from which additional contraints 
may ultimately be derived ||. The present study was 



motivated by an isospin dependent Boltzmann-Uchling- 
Uhlenbeck (BUU) transport model calculations of Ref. 8 
for 112 Sn+ 112 Sn and 124 Sn+ 124 Sn collisions. Using pro- 
ton and neutron densities calculated from the non-linear 
relativistic mean-field theories as inputs, the dynamics of 
nucleon-nucleon collisions are calcuated. These calcula- 
tions utilize nucleonic mean field potentials correspond- 
ing to an equation-of-state that can be expressed (at zero 
temperature) in terms of the mean energy of a nucleon 
in nuclear matter consisting of a sum of terms from a 
symmetric part and an asymmetric part: 



E(p,5)=E(p,S = 0) + S(p) 6 2 



(1) 



where the mean-field component of the symmetry energy 
is given as a product of the square of the local asymmetry 
and a density dependent factor S(p), which is sometimes 
described at sub-saturation density by a power law S(p) 
= Cu 7 where u= The uncertainty in 7 (0.5< 7 <1.6) 
reflects uncertainties in the nucleon effective interaction 
that require better experimental constraints [fj. Possi- 
ble density dependences of the asymmetry term of the 
EOS - defined to be the relationship between the density 
and the energy (or pressure, P=p 2 d(E/A)/dp at constant 
entropy) range from "asy-stiff" expressions in which the 
energy and pressure increases more sharply as a func- 
tion of density, "asy-soft" in which the pressure does not 
increase as sharply. 

Calculations of neutron and proton spectra, for mo- 
mentum independent mean field potentials with 7=0.5 
and 7 wl.6 reveal a strong sensitivity of the neutron- 
proton spectral ratio R„/ p (E cm ) = (dY„(E cm )/dE cm ) 
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/ (dY p (E cm )/dE cm ) of transversely emitted nucleons to 
the density dependence of the asymmetry energy, where 
E cm is the particle energy in the center-of-mass frame. 
For a strongly density-dependent - or "asy-stiff" - sys- 
tem (with 7 «1.6) at sub-saturation densities, the ratio 
changes little with system asymmetry S due to the weak- 
ness of the potential at u<l compared to the Coulomb 
interaction. However, for a weakly density-dependent - 
or "asy-soft" - system (with 7=0.5), R„/ p (E cm ) depends 
strongly on S, and a larger emission ratio is found for the 
124 Sn+ 124 Sn reaction than for the 112 Sn+ 112 Sn reaction, 
reflecting the preferential expulsion of neutrons from a 
neutron rich system due to the symmetry energy |§j]. 

Beams of 112 Sn and 124 Sn at E/A=50MeV were pro- 
duced in the K1200 cyclotron at the National Super- 
conducting Cyclotron Laboratory (NSCL). These beams 
were incident on 112 Sn and 124 Sn targets located within 
a thin-walled aluminum target chamber. The chamber 
contained a central charged-particle multiplicity detec- 
tor, a fast plastic scintillator array at forward angles, 
and three charged-particle telescopes. Two large area 
neutron walls were place outside the chamber. Each wall 
had active areas of about 2 m square with each consist- 
ing of 25 horizontally stacked 2 m long quartz glass tubes 
filled with NE213 liquid scintillator 10] . One wall was 
placed at a distance of 4.59 m covering forward polar an- 
gles of 8° < 9iab <34°, and the other wall was placed 
at a distance of 5.79 m, covering forward polar angles of 
42° < 9iab <62°. Associated pulse shape discrimination 
techniques 0] were used to distinguish neutrons from 7- 
rays. During the experiment, shadow bars were inserted 
between the scattering chamber and the neutron walls to 
assess the background yield of neutrons scattered from 
the floor and walls of the experimental vault. Neutron 
and proton emission was observed at angles between 70° 
and 110° in the CM to suppress the contribution from 
decays of the projectile- like fragment (PLF). 

Neutron energies were measured by time of flight rela- 
tive to the trigger from a segmented BC404 start scintil- 
lator array that was placed inside the chamber at a dis- 
tance of 10 cm downstream from the target. This array 
consists of four 3 mm thick plastic paddles, each trape- 
zoidal in shape, arranged at forward angles covering an 
angular region of 5° to 60°. Each paddle in the array pro- 
vided an azimuthal angular coverage of 72°. A "missing" 
fifth paddle provided a clear flight path from the target 
to the charged-particle detector array described below. 
The time resolution of the plastic scintillator array was 
less than 300 ps FWHM and the neutron array time res- 
olution was measured to be about 700 ps FWHM. 

Isotopically resolved charged particles with 1<Z<3 
were measured using three telescopes of the Large- Area 
Silicon-Strip Array (LASSA) [H|. Each telescope con- 
sists of a 500 ^m double-sided silicon strip detector 
(DSSD) followed by four 6 cm thick CsI(Tl) scintilla- 
tor detectors arranged in quadrants. Each side of the 



DSSD consists of 16 strips per side with a strip pitch 
of 3 mm. The orthogonal strips on the front and back 
sides of each strip detector subdivided the 25 cm 2 square 
surface area into a Cartesian grid of 256 pixels. The 
LASSA telescopes were arranged in the horizontal plane 
opposite the neutron array at a distance of 20 cm from 
the target for a total polar angle coverage ranging from 
15° < 6 lab <60°. These telescopes were calibrated with 
an accuracy of about 4% using recoil protons scattered 
from a CH 2 target by an 18 beam with E/A=50 MeV. 

Central collisions were isolated using overall multiplic- 
ity as an indicator of the impact parameter fl2T ]. The 
impact parameter was measured using the Washington 
University MicroBall 0|. In its standard configuration, 
the MicroBall consists of 95 closely packed CsI(Tl) scin- 
tillators arranged in 9 rings, covering a polar angle of 4° 
to 172° providing about 97% of 4tt coverage. Scintilla- 
tion light is collected by silicon photodiodes, and particle 
identification for 1,2:3 H, 3,4 He, Li, Be, and B is accom- 
plished using associated pulse-shape discrimination elec- 
tronics. For this particular experiment, charged particles 
emitted at polar angles of 60° < 9 <172° were detected 
in 55 elements of the partial MicroBall array. Assuming 
the charged particle multiplicity decreases monotonically 
with impact parameter, a "reduced" im pac t parameter 
b=b/b max was determined following Ref. Il2t here, events 
with a charged-particle multiplicity of Nc >7 corre- 
sponding to b <0.2 were selected. By comparison to sim- 
ulations of the same collisions with the Anti- Symmetrized 
Molecular Dynamics (AMD) computer code we as- 
sessed that this data set corresponded to events uni- 
formly distributed over impact parameters of b<5 fm. 
Comparison to the AMD is made becuase the BUU code 
does not accurately accomodate coalescence, and hence, 
the accuracy of the charged-particle multiplicity at lower 
energies suffers. 

Neutron-proton transverse emission ratios in the CM 
were compared to predictions given in Ref. 8 for the 
two systems - 124 Sn+ 124 Sn (5=0.19) and 112 Sn+ 112 Sn 
(5=0.11). These ratios are predicted to vary with 8 ac- 
cording to the stiffness of the EOS. For a relatively asy- 
soft EOS, proton emission is supressed relative to the 
Coulomb repulsion at the sub-nuclear densities studied, 
while proton emission becomes comparable to neutron 
emission for a relatively asy-stiff EOS. A single measure 
of the stiffness of the asymmetry term is given by the 
double ratio 1R124/1R112 = (dn n /dn p )i24/(dn„/dn p )ii2 
of A=l products for the two reactions indicated by the 
mass of the projectile (and target) in the subscript [l5| . 
The double ratio is useful in that - to first order - detector 
efficiencies, Coulomb effects, and systematic errors cancel 
out, increasing the sensitivity over that of the single ratio 
to the symmetry energy. Ratios of the nuclei 3 H- 3 He are 
also studied since these are predicted to be sensitive to 
the density dependence of the nuclear asymmetry term 
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FIG. 1: Ratios of t/ 3 He and n/p for the 112 Sn+ 112 Sn reaction 
(a) and the 124 Sn+ 124 Sn reaction (b) as a function of energy 
per emitted nucleon for particles emitted between 70° and 
110° in the center-of-mass. In the top plot, error bars are 
smaller than the points. 




40 50 

E=.„/A(MeV) 



FIG. 2: Double ratios 1R124/1R112 for free nucleons emitted 
in each reaction compared to the double ratios of free nucleon 
yields calculated from the BUU transport calculations of Ref. 
iSi. The filled circles correspond to double ratios of yields of 
tranversely emitted free nucleons and the open squares corre- 
spond to all nucleons including those bound in clusters. 



|l(j| . Neutron and charged particle energy spectra were 
obtained for energies up to about 100 MeV in the CM. 
Neutron-proton and 3 H- 3 He transverse emission ratios 
are shown in Figure Q] for the two reactions studied. It 
is noted that the ratios corresponding to the A=3 mir- 
ror nuclei are very similar to those of free nucleons above 
20 MeV/ A with the most discrepancy at low energies 
where binding energy and Coulomb effects become con- 
siderable. The A=3 mirror nuclei ratios are within 10% 
of those obtained in recent experiments performed by Liu 
et al. using the same reactions at the same energy 0. 
Both ratios are lower in the 112 Sn+ 112 Sn reaction than in 



the 124 Sn+ 124 Sn reaction. Qualitatively, the results pre- 
sented in Figure ^ do not necessarily indicate an asy-soft 
EOS, as one naturally expects a higher neutron-proton 
ratio from the 124 Sn+ 124 Sn system. However, given the 
similarity of the ratios at higher center-of-mass energies, 
this provides evidence that mirror nuclei transverse emis- 
sion ratios may provide a good measure of the density 
dependence of the asymmetry energy. 

Figurenpresents results from an analysis in which clus- 
ter emission and free nucleon emission are treated inde- 
pendently. However, cluster production may be less sen- 
sitive to the asymmetry potential for higher energy clus- 
ters ^||. More importantly for this analysis, the effects 
of evaporation from light clusters and the PLF as well 
as effects from coalescence may not be well understood 
in the transport calculation. In order to explore the ef- 
fects of cluster yields from fragment decays, a coalescent- 
invariant analysis has been conducted in which isospin 
emission ratios for a reaction f2 are taken as the ratios of 
all emitted nucleons - both free and bound: 



p _Y n _ ^ i N i Y i 



(2) 



where the summation is over all species i with charge 
Zj, neutron number Nj, and yield Y, up to mass A=16 
and charge Z=6. Because the current experiment only 
measured complete energy spectra over the entire range 
of particle energies for nuclei with Z<2, the fragment 
yields for this analysis were taken from the experiment 
of Liu et al. in which the identical reactions were studied 
and fragment spectra for Z>2 were measured 0. 

In order to compare the nucleon transverse emission 
ratios from each reaction, the double ratio was used [lij . 
The double ratios for free neutrons and protons are shown 
in Figure Also shown are the double ratios calculated 
from a BUU transport calculation for two asymmetry re- 
lationships. In this calculation, a nucleon is considered 
free if it is not correlated with another nucleon within a 
spatial distance and a momentum distance of 3 fm and 
300 MeV/c respectively at a time of 200 fm/c after the 
initial contact of the two reacting nuclei . The BUU re- 
sults correspond to an asy-stiff relationship, with a bulk 
symmetry energy defined by S'cx2u 2 /(w-(-l)~u 16 , (shown 
in the shaded region) and an asy-soft relationship with 
a bulk symmetry energy Socu 1 / 2 (shown in the cross- 
hatched region). The BUU results shown are for central 
collisions with b<5 fm, which is identical to the impact 
parameter for the shown experimental results. Relevant 
nucleon emission is only displayed for energies above 20 
MeV/ A as emission at lower energies is dominated by 
cluster emission and is not well modeled by BUU calcula- 
tions. Coulomb effects at low energies can also adversely 
affect the comparison. 

Care was taken to eliminate some of the experimen- 
tal uncertainties, as well as those associated with the 
BUU calculation. First, because BUU calculations do 
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not include evaporation (while the experiment necessar- 
ily does), nucleon emission from the PLF is not included 
in the theoretical results. This is overcome by examin- 
ing transverse nucleon emission, because PLF emission is 
primarily at forward angles. Second, effects from coales- 
cence are examined via the coalescent invariant analysis 
of Equation [21 in which all nuclcons - free and bound - 
are counted in the ratios. The results of this analysis are 
shown in Figure [3 Because the results of the coalescent- 
invariant analysis closely match those of free nucleons at 
energies above 25 MeV, it is concluded that the trans- 
verse emission at higher energies is entirely free nucleon 
emission, and effects from coalescence become negligible. 
The only remaining uncertainty is due to the fact that no 
angular cuts are made in the BUU calculation due to lim- 
ited statistics. Currently, BUU calculations are planned 
to improve the statistics near 90° in the center-of-mass 
so that such angular cuts will be possible. These will be 
available in the future, and preliminary results indicate 
that the inclusion of angular cuts has little effect on the 
overall results [T?) ]. 

The experimentally determined double ratio (the ra- 
tio of Y„/Yp for the 124 Sn+ 124 Sn reaction to that for 
the 112 Sn+ 112 Sn reaction), which may be understood as 
the change in the proton repulsion (relative to the neu- 
tron repulsion) with system asymmetry, appears to cor- 
respond more closely to an asymmetry potential in which 
u < F(u) < it 1 / 2 . In these results, the relative proton 
emission with respect to to the neutron emission increases 
by over a factor of two as the system asymmetry goes 
from 0.11 ( 112 Sn+ 112 Sn) to 0.19 ( 124 Sn+ 124 Sn). This 
can be understood by considering the proton and neu- 
tron mean field potentials corresponding to strongly and 
weakly density dependent asymmetry energies at sub- 
saturation densities. For the densities less than po stud- 
ied here the proton (neutron) attractive (repulsive) po- 
tential for an asy-soft EOS increases more sharply with 
density than for an asy-stiff EOS. The net result is that 
the proton attraction competes more with the Coulomb 
repulsion in the asy-soft case as the system asymmetry 
changes, producing the larger double ratio shown in Fig- 
ure El due to the repressed proton emission and the en- 
hanced neutron emission, compared to the asy-stiff case 
in which the proton emission changes little as the asym- 
metry changes. These results are consistent with pre- 
dictions of the 208 Pb skin thickness and with predictions 
from isospin diffusion, in which the deduced asymmetry 
energy S is roughly proportional to u om \\A Il9i|. 

Nucleon transverse emission has been studied as a sen- 
sitive observable of the asymmetry term of the nuclear 
EOS. Comparison to BUU transport calculations indi- 
cate a somewhat asy-soft EOS. Measurements of total 
nucleon emission ratios - including those bound in clus- 
ters - is suggested as a way of compensating for possible 
ambiguities in coalescence calculations, which are largest 
at the lowest energies. The coalescent-invariant results 



presented can be compared to transport calculations as 
confirmation of current theoretical understanding of the 
density dependence of the nuclear asymmetry term. 

Further experimentation should also concentrate on 
observables which are predicted to be sensitive to the 
high-density behavior of the asymmetry term. This in- 
cludes not only transverse nucleon emission as described 
in the current research, but may also include total pion 
production yields . Future neutron- proton transverse 
emission is planned at the NSCL for slightly higher den- 
sities, while facilities like the RI Beam Factory may be 
suitable for explorations at larger 8 and at densities sev- 
eral times Pq. 
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